Dark States and Interferences in Cascade Transitions of Ultra-Cold Atoms in a Cavity 
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We examine the competition among one- and two-photon processes in an ultra-cold, three-level 
atom undergoing cascade transitions as a result of its interaction with a bimodal cavity. We show 
parameter domains where two-photon transitions are dominant and also study the effect of two- 
photon emission on the mazer action in the cavity. The two-photon emission leads to the loss of 
detailed balance and therefore we obtain the photon statistics of the cavity field by the numerical 
integration of the master equation. The photon distribution in each cavity mode exhibits sub- 
and super- Poissonian behaviors depending on the strength of atom-field coupling. The photon 
distribution becomes identical to a Poisson distribution when the atom-field coupling strengths of 
the modes are equal. 

PACS number(s): 42.50.Vk, 42.50.Dv, 03.75.-b 
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I. INTRODUCTION 



The dynamics of atom-field coupling in high quality cavities has been an interesting tool to verify the predictions 
of quantum electrodynamics (QED) on the radiation-matter interaction [0-^. The properties of atoms in cavities can 
be modified in a controlled way by the proper design of cavity geometry and its quality factors. For example, the 
spontaneous emission of an atom can be enhanced or inhibited inside a cavity Cavities with high quality factors 
are especially attractive for studying the regime of strong interaction between the atoms and the quantized fields. 
Interesting effects of strong interaction such as collapse and revival of Rabi oscillations, vacuum Rabi splitting, atom- 
atom and atom-field entanglements have been predicted and observed (^jH) . The generation of atom-field entangled 
states in cavity QED finds many applications in quantum information processing, logic gates etc. A very important 
result of strong interaction was the micromaser where the excited atoms go successively through a high quality cavity 
in a time much smaller than the characteristic decay times. This leads to build up of a steady state field in the 
cavity which has properties quite different from laser fields. The theoretical works on micromaser suggested 
the sub-Poissonian photon statistics which is purely due to the quantum nature of the field. The operation of the 
micromaser has been realized and Weidinger et al verified the existence of trapping states in the micromaser . 
The two-photon micromaser as well as the microlaser were also realized pO| , pl| . The theory of two-mode micromaser 
operating on three-level systems has been developed [|l2| . 

In all these works, the external motion of the atom is treated classically while its interaction with the field is 
quantum mechanical. By considering ultra-cold atoms, Scully et al discovered that the quantum treatment of 
center-of-mass (cm.) motion of an atom leads to a completely new kind of induced emission in a cavity. The dressed 
state analysis shows that the cavity field can act like a quantum mechanical potential for an ultra-cold, incident atom. 
The nature of the potentials depends on the mode profile of the cavity as well as the atom-field coupling strengths 
p^ , p^ . Thus the atom besides changing its electronic states can be either reflected or transmitted through the cavity. 
This reflection or transmission of the atoms is very similar to those of a particle interacting with potential barriers 
or wells. We have recently shown that tunneling of an ultra-cold, two-level atom in the excited state through two 
successive cavities exhibits transmission resonances similar to those of an electron tunneling through semiconductor 
double barriers [|l6|. The reflection and transmission through cavity induced potentials give rise to many interesting 
features in the photon statistics and the spectrum of the micromaser pumped by ultra-cold atoms Further, 
the works on quantum treatment of the atomic motion have been extended to study the interaction of two- and 
three- level atoms with a single mode fleld Considerin g th e quantized motion of A-type three- level atoms, the 
photon statistics of a two- mode micromaser was discussed [|l9| . In the A-type scheme, an excited atom can make 
only one photon transition to either of the two ground levels. In this paper, we examine the photon statistics of a 
two-mode micromaser pumped by ultra-cold, cascade three-level atoms. The excited atom in the cascade transition 
emits photons sequentially in each of the two cavity modes. Thus, the field mode which is resonant with the upper 
transition of the S-type atoms, gets populated by both one and two-photon transitions from the excited state. We 
discuss the effect of two-photon emissions from the excited atoms on the photon statistics of the micromaser. We 
make use of methods similar to those in Ref . . 

The organization of the paper is as follows. In Sec. II, we describe the interaction of ultra-cold, S-type three-level 
atoms with the bimodal field in the cavity. We show that the field can always induce two-photon transition in an 
excited atom when the corresponding one-photon transition is forbidden. In Sec. Ill, we study the mazer action and 
derive the master equation for the reduced density operator of the field in the cavity. In Sec. IV, the numerical results 
are presented for the steady state photon distribution in each mode of the bimodal field. We find that the photon 
distribution in each mode exhibits Poissonian statistics if the atom-field coupling strengths of the cavity modes are 
equal and when the two-photon transitions from the excited atoms are dominant. For different coupling strengths of 
the cavity modes, the dominant two-photon effects in transitions can also lead to sub- and super-Poissonian statistics 
of photons in the steady state field. Finally, we consider the mazer action when the pumping atoms have kinetic 
energies close to the vacuum coupling energy. In this case, one-photon effects also become important and these lead 
to different photon statistics for the two cavity modes. 



II. ONE AND TWO PHOTON PROCESSES IN ULTRACOLD ATOMS IN A CAVITY 



We consider a bimodal cavity of length L pumped steadily by a beam of ultra-cold three-level atoms in the cascade 
configuration. The scheme of our model is shown in Fig. 1. The transitions a —^ bi and bi —^ 62 are dipole allowed 
while the direct transition a — > 62 is dipole forbidden. Thus the atom in the excited level a can reach the ground level 
62 only through the two-photon transition a bi — > 62 • The frequencies of the two cavity modes 1 and 2 are tuned 
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to those of atomic transitions a ^ bi and bi —^ 62 respectively. The Hamiltonian describing this resonant atom-field 
interaction including the quantized motion of center-of-mass of the atoms along the z direction is given by 

H ^ Ha + Hf + Haf , (1) 

where Ha{Hf) is the Hamiltonian of the free atom (field) and Haf is the interaction Hamiltonian describing the 
atom-field interaction in the dipole and the rotating wave approximations: 

2 2 
Ha = ^ + hna\a){a\ + Y,^^bJbo. )(6a| , 

a— 1 



Hf = 



E 



huJatt'ltta , 



Haf = ngi{ai \a){bi\ + \bi){a\ al) + ^2(02 |&i)(^'2| + |^'2)(6i| 4) 



(2) 



The operator — 0,61,62) gives the projection on to the state |j) with energy hflj. The operators |a)(6i| 

and |6i)(62| describe the atomic transitions from the upper and lower levels to the middle level. The operators 
(a]j) annihilate (create) a photon in the modes a — 1,2 with resonance frequencies uji = fla — and UJ2 = — 
respectively. The first and second terms in the interaction operator Haf represents the action of fields 1 and 2 
of the cavity on the upper (a <^ 61) and the lower (61 62) transitions respectively. The parameters ga are the 
corresponding atom-field coupling constants and m is the atomic mass. The parameters ga are dependent on z through 
the mode function of the cavity. 

In the interaction picture, the Hamiltonian ([l]) of the atom-field system reads 



Hr = 



2m 



H 



AF ■ 



(3) 



It is useful to expand the interaction Hamiltonia n Haf hr its diagonal basis. The operator Haf has eigenstates 

{ni 



/>°,+i,„,+i), \<Pt,+i,n2+i) witl^ eigenvalues 0, ±h^y g1{z)(ni + 1) + g'i(z){n2 + I), respectively, where 



+ I,n2 + 1/ 



V5l("l + l)+52("2 + l) 



|a, 711,712) 



1 



V5iK + f)+ff|(ri2 + f) 



162,711 -f I, 712 4- 1) 



l</'rti-|-l,n2 + l) — ^ 



|a, Til, 112) ± |6i,7ii -I- 1,712) 



v/ff?(ni + l)+.gi("2 + f) 

:|62,71l + 1,712 + 1) 



g2V"-2 + 1 



y/gfini + l)+g^{n2 + l) 



(4) 



The interaction operator Haf and its eigenstates |0°^+i „2+i), l^^^+i „2-i-i) depend on the position z through the 
coupling strengths 51(2), 32(2). Thus, it is in general difficult to carry out the time evolution of an atom-field state 
governed by the Hamiltonian (^) for the quantized motion of atoms. So, for simplicity, we work with the mesa 
mode functions ga{z) = ga 9{z)9(L — z) which represent the constant field modes in the cavity. In this case, the 
eigenstates of the interaction are independent of atomic position inside the cavity and the atomic motion sees free 
particle evolution in the dark state of interaction |0°j_|_i „2+i)- The effect of atom's interaction with the cavity on 
its external motion can be realized only in the dressed state n^+i) components of the initial atom-field state. 

We need to consider the initial atom-field state to be |a, 711,712), i.e., the atom is in the excited state and the cavity 
contains (711,712) photons in the modes (1,2) initially. This initial state can be expanded in the dressed state basis as 



a,7ll,7l2> 



gi^{ni + l)/2 



^/gi{ni + l)+gi{n2 + l) 



+ 



.92\/"2 + 1 



V5?K + l)+5i("2 + l 



(I<^ni + l,n2+l) + l'/'ni+l,n2 + l)) 



■I'^ni+l,n2 + l) 



(5) 
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The time evolution of this initial state can be found by expanding the combined state of atom-cavity system as 

l*(2,i)> = X+(^>*)l0ni + l,„2 + l> +X-(^,OI'/'r^i+l,„2 + l) +Xol'?^"l+l,«2 + l) ' (6) 

then the time dependent Schrodinger equation becomes 

dXa{z,t) , / ,^ , r, 

iTl — = haXa{z,t) , q; = ±,0. (7) 

Ot 

Here, h± = p1/2m ± hy/ gl{ni + 1) + g|(»T.2 + 1), ho = pj./2m are operators acting in the space of center of mass 
variables. Thus, the effect of the cavity with fixed number of photons produces potential terms in corresponding 
to the dressed states If/'^^^.i „2+i) as discussed in Rcf. The barrier and well potentials induced by the interaction 
for the atomic motion in the states |(/)^^_|_]^ n2+i) then displayed as in Fig. 2. It is also important to note that the 
external motion of atom experiences free time evolution in the dark state If/i^j+i „2+i) for the mesa mode distribution 
of the cavity fields. Denoting the reflection and transmission amplitudes as p^^ t^^ „^ for the potential barrier- well 
problem of the dressed states \4'ti+i.n2+i) ^ respectively, we have 

Pni,n2 = *^ni,n2 sin(fc,ti ,„2 i) exp(ifcL)T^^ „^ , (8) 

^«i,"2 = exp(-iA:L)[cos(fc^^ „^L) - i^t^,n2 sH^t^^n^L)]'^ , (9) 



1 /fc 



2\ k fc± , 



ktn. = y (fc' T ^v'5?("l+l)+.9iK + l)) , (11) 

where hk is the atomic cm. momentum and L is the length of the cavity. It is to be noted that the strengths of 
the barrier - well potentials (potential energy term in k^^ depend on the coupling constants gi, g2 as well as the 
occupation numbers ni, n2 of the photons in the cavity. 

We consider the initial wave packet of the moving free atom to be ip{z,t) — exp (^~iplt/2mlij J dkA{k)e^''^ ~ 

J dkA{k)e^^^^'' /2m)t^ikz^ rp-j^g Fourier amplitudes A{k) are adjusted such that the peak of the incident wave packet 
|V'(0, at the entry (z — 0) of the cavity occurs at time t = 0. The combined state of the atom-cavity system at 
the initial time i = is therefore, 

|*(z,0)) =^(z,0)|a,ni,n2) . (12) 

The wave function of the atom-field system at time t is found by solving the Eqs. (||) and subject to the above 
initial condition : 

\^iz,t)) = J dkA{k)e-<''''"/'"'y {[Ra,n^^nMe-''''H-z) +Ta,n,,nMe^^^ \a,n,,n2) 

+ [Rb,,m+l,nMe-"''0i-z)+n,^n, + l,nAky''"O{z - L)] \b,,n, + 1,712) 

+ [i?b2,«l + l,n2 + l(fc)e~''''^(-2)+T^fc2,ni + l,„2-hl(fc)e*'''^?(2-i)] 162,^1 + 1,712 + 1)} , (13) 



where 



p _ g?("l + l) ( + , - N 

^a,m,u2 - 2{gl(ni + l)+ gl{n2 + 1)) ^^"^'"^ ^"^'"^^ ' 

'"'"-"^ 2 {glim + 1) + gl[n2 + 1)) ^ ^ -'-^ ^ {gl{n, + 1) + gl{n2 + 1)) ' ^ ^ 
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are the reflection and transmission amplitudes for the excited state of the atom with the initial (711,77,2) photons 
remaining in the two cavity modes and 



P giV("i + 1) f + _ - N 

2v(5rK + 1) + 52(^2 + 1)) 

_ giVini + 1) . + _ _ X , . 

2v(.9r("i + 1) + 52(^^2 + 1)) 

are the probability amplitudes that the excited atom goes to the state \bi) and emits a photon in mode 1 while getting 
reflected and transmitted respectively. Similarly, the excited atom is reflected or transmitted and emits a photon in 
both the cavity modes while making a transition to the ground state I62) via the middle state with probability 
amplitudes 

p _ .9152^(771 + l){n2 + 1) . + N 

nb2,n, + Un2 + l - 2 (^2(„^ + 1) + ^2(„2 + 1)) {Pm.n^ + Pm.n^) > 

rj. gig2v/(7ti + i)(7i2 +1) _ . gig2\/{ni + i)in2 + 1) 

2(.gf(7T.i + 1) +g^(?i2 + 1)) ' + 1) +5^(772 + 1)) 

It is clear from the above equations that the effects of the barrier and well potentials induced by the dressed states 
I'/'ni+i n2+i) ^^'^ coherently in either the reflection or transmission of the atom. The additive term to the barrier-well 
amplitudes in Eqs. and ( [l6| ) comes from the contribution of the dark state „2+i) in the initial state 

expansion Eq. (|^). Since the dark state is orthogonal to the middle state it influences only the two-photon 
emissions and not the one-photon emissions of the excited atom. An important feature here is that the two-photon 
transition can always be induced by the field when the one-photon transition is forbidden. This can be seen by 
examining the probabilities for different states of the atom. When an initially excited atom is incident upon the 
cavity containing (71,1,712) photons in the two cavity modes (1,2), respectively, then from Eqs. (|l5|) and ( p^ ) the 
probability that the atom makes a one-photon transition to the state |6i) with the emission of a photon in mode 1 is 

^'ni, 712(0 ^1) = |-Rbi,ni-|-l,n2 1^ + l^bi ,ni + I,n2 I ^ J (17) 

and the probability that the atom makes a two-photon transition to the state I&2) with the emission of a photon in 
each of the modes 1 and 2 of the cavity is 

Pni,n2{0' ^ ^2) = |-Rb2,ni + l,n2 + l|^ + \Tb2Mi + l,n2 + l\'^ ■ (18) 

When (72 = 0, the lower transition bi — s- &2 is forbidden and hence the two-photon transition probability Pni,n2 (^ ~* ^2) 
vanishes. Then, the upper transition a — > 61 behaves like a two-level atom interacting with the mode 1 of the cavity 
field. In Fig. 3, we compare the photon emission probabilities of excited, two-level (52 = 0) and three-level (52 7^ 0) 
atoms when the cavity is initially in vacuum state. We scale the parameters in terms of a wave number k which 
is defined by the vacuum coupling energy hgi = h^n^ /2m of the two-level atom. Note that the coupling strength 
gi = 2tt X 10 MHz of a *^Rb atom corresponds to 155 /im length of the cavity for the parameter kL — 200007r. The 
temperature of the atom is of the order of 10~^ K for the mean momentum k/K = 0.01. For ultra-cold, incident 
atoms [k/K « 1), the graph in Fig. 3(a) shows that the one-photon emission probability exhibits resonances as a 
function of length of the cavity for both two-level and three-level atoms. In the three-level atom, the resonances of the 
one-photon emission probability occur at values of K.L different from those of two-level atom. The two-photon emission 
probability shows maxima and minima at the resonance positions of the one-photon transition. This behavior arises 
from the interference of the contributions coming from different dressed states in the initial state expansion Eq. (|^). 
The transmission probability \Tb^^rii+i,n2+if' obtained from Eq. ( p^ ) has an interference term proportional to the 
phase of the amplitude (t,+ „^ -I- t,~ „^ ) . This term can be constructive or destructive which leads to the enhancement 
or reduction of the two-photon emission probability. From the graph, we also see that Po,o(i — * ^2) ~ 0.32 when 
Pq,o{0' — * ^1) ~ 0. This implies that the probability of two-photon emission is not the product of the probabilities for 
single photon emission. This also suggests that both the field modes of the bimodal cavity can be amplified together 
through the two-photon transition of an excited atom. This feature is absent in the case of two-mode micromaser 
pumped by ultra-cold, A-type three-level atoms ]l9| . In the A-scheme micromazer, the excited atom can make only 
one-photon transition to either of the two ground levels. Therefore, both the cavity modes can not be populated 
sequentially by the photon emission from the same atom. Next, we compare our results for ultra-cold atoms with 
that of fast, incident atoms in the 2- type configuration. In the case of fast, incident atoms [k/n » 1), both onc- 
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and two- photon emission probabilities exhibit Rabi oscillations as a function of the length of the cavity as seen from 
Fig. 3(b). The Rabi frequency of oscillation for one-photon emission is twice that of the two-photon emission. These 
features resemble exactly the results in the usual Jaynes-Cummings (JC) model where one neglects the quantization 
of the atomic motion in atom-cavity interaction Thus, neglecting the kinetic energy operator in the Hamiltonian 
the time evolution of the initial atom-field state |a, rti, 712) gives for the photon emission probabilities in the JC 
model: 

P„..n.(a-.6i) = ^i^^^i±^sin^(f^r) , 



Pn,,n2{a 02) = ^ sm (S2r/2) , 

where il = glini + 1) + .91(^^2 + 1) is the Rabi frequency and r is the interaction time of the atom with the cavity. 
Note that when Qt is an odd number multiple of tt, the one-photon transition is forbidden and the two-photon 
transition probability becomes maximum consistent with the results for fast atoms in Fig. 3. 



III. MAZER ACTION IN A BIMODAL CAVITY - BASIC MASTER EQUATION FOR THE CAVITY 

FIELD 

In this section, we consider the mazer action of ultracold atoms in a bimodal cavity and derive the master equation 
for the cavity field assuming that excited atoms are pumped steadily into the cavity. We model the random pumping 
of the atoms into the cavity by a Poissonian process with an average rate of pumping r. The fiux of the incident 
atoms is adjusted so that only one atom interacts with the cavity field at a time. We neglect the cavity field damping 
during the time an atom interacts with the field. Since the field in the cavity changes with the passage of each atom, 
we need to know the time evolution of the atom-field state for a general initial state of the cavity field. The wave 
function of the initial atom-field system is now given by 

^nja, 711,712). (19) 

ni,n2 

Carrying out the time evolution for this initial state using Eqs. (||)-(0), the state of atom-field system after the 
interaction is given by 

I ^{z,t)) = J dfcA(fc)e~'('''='/2")* 

00 

J2 [T^a.nun2{k)e~"'""0{-z) | fl, Tii , Us) + Ta,„i (fc)e*'=^ 6'(z - L) | 0,711,712) 

ni ,n2— 

+ 7^b„„,+l,„,(fc)e-*'="^?(-z) I bi,ni + l,n2) +Tb,^n,+l.^Aky'''eiz - L) (20) 
I bi,ni + 1,712) +7^^,2^„l+l,„2+l(fc)e"*''^6'(-z) | 62,111 + 1,713 + 1) + Tt2,ni+i,n2+i{k) 
e''''-e{z~L) I 62,711 + 1,112 + 1)] , 

where 

'T^a.ni ,n2 (^) ,n2 -^a,ni ,n2 (^) : 

^^,ni,n2(^) — ,n2 '^a,ni .712 (^) ' (^-i) 

are the probability amplitudes for refiection (or) transmission of the atom in the upper state \a) with (rii, 712) photons 
in the cavity field and similarly, the atom is refiected (or) transmitted when the atom-field state is |6i, ni + 1, 712) (or) 
1 62, Til + 1, 112 + 1) with amplitudes 

"^bi ,ni + 1,712 (^) — ,712 -^61 ,ni + 1,712 (^) ' 

^1 ,711 + 1,772 (^) ^Tli ,712 -^61 ,711 + 1,712 (^) ' 

'^f)2 ,711 + 1,772 + 1 (fc) — C'711, 712^62, 711 + 1, 7l2 + l(fc) ) (22) 
^2,7ll + l,7l2 + l(^) — ^^711, 71221)2, "1 + 1, 7l2 + l(^) ■ 
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The time evolution of the reduced density operator of the field in the interaction picture is then given in the coarse 
graining method [Q to be 

p{t) = r5p[t) + Lp{t) , (23) 

where 5p{t) = p{t) — p{<S) is the change in the reduced density operator of the field due to the passage of a single 
atom in the excited state. This can be obtained by forming the atom-field density matrix using Eqs. (|l|) - (H) and 
tracing over external and internal degrees of freedom of the atom. Field damping and the effect of thermal photons 
are described by the Liouville operator 

Lp = ^Ci(nfcj + l){2aipa\ - a\aip - pa\ai) 
+ ^CiUb^ {2a\pai - aia\p - paia\) 

+ -jC2{nb^ + l){2a2pa\ - a\a2P ~ pa\"''2) (24) 
+ ^^271^2(20^^02 - a2a\p- pa2a\) . 

Here Tif,^ is the number of thermal photons in mode a and Ca is the damping rate of this mode. Using Eqs. ( p3| ) and 
( p^ wc obtain the equation governing the time evolution of density matrix elements, 

p{ni,n2]n\,n'2) ^ r {{Ra,nun^K,n^' ,n^' +Ta,nunJ'a,m',n^' " l)p(ril , ; ^I'l , 4 ) 

+ (^bi,ni,n2^bi,ni',n2' + ^bl ,"1 ,"2 ^b* ,ni ' ,712' ~ 1; '^2] - 1, n2) 

+ (^f)2,ni,«2^b2,«i',n2' + '^&2.ni,n2^b*2,«l',«2')/'('^l ^ 1' ^"^ ^ ^5 "l' ~ 1' "■2' " 1)} 

+ ^C-iKi + l)[2y^ (rii + l)(n; + l)p(?ii + 1, ^2; < + 1, n'2) 

- (ni + n[)p{ni,n2; n[,n'2)] 

+ ^Cinfei [2^ nin'^p{ni - l,n2;n[ - 1,122) (25) 

- (ni + n[ + 2)p{ni,n2;n[,n2)] 

+ \c2{nb., + I)[2^{n2 + l){n'2 + l)p{ni,n2 + l;n'i,n'2 + 1) 

- (^2 + ri2)p(ni, ?i2; n'^, ^3)] 



+ -(72^62 [2 V n2n'2p{ni,n2 ~ l;n[,n2 - 1) 
- (n2 + 7i2 + 2)p{ni,n2; 713)] . 

The diagonal elements of the density matrix P(rti,rt2) = p(rti, ^2; ni, 712) which gives the joint distribution of 
photons in the two cavity modes, obeys the following equation: 

P(ni, ^2) = -G'b,,„i,„2P(ni, 712) + Gbi,„i-i,«2P(n.i - 1, ^2) 

- Gb2,ntjn2P{ni,n2) + Gb2,ni-1,„2-1-P("-1 - 1, 7l2 - 1) 

+ Ci(nb, +l)[(ni + l)P(ni + l,n2)-niP(ni,n2)] (26) 

+ CiUb^ [niP{ni - l,n2) - (rii + l)P(ni,n2)] 

+ C2(nb2 + 1) [("2 + l)P(ni, ^2 + 1) - n2F(ni, ^2)] 

+ C2rib2 [n2P{ni,n2 - 1) - (7^2 + l)P(7ii, 712)] , 

where Gbi.ni,n2 — i'Pni,n2{'^ ^ ^1) E^J^d Gb2,ni.n2 — fPni.7i2{^ ^2) are the gain coefficients for the atomic transitions 
with Pni,n2{o- — bi) and P„j^„2(a — > ^2) as defined in Eqs. ( |l7| ) and (|l8|). This is the master equation for the two- 
mode micromaser describing the time evolution of photon distribution in the cavity. This equation behaves similar to 
a rate equation for the probability and a simple physical meaning can be given to each term on the right hand side in 
terms of inflow and outflow of probabilities. The first and second terms in the equation gives the effect of one-photon 
transitions while the third and fourth terms correspond to the two-photon transitions of the excited atoms. 
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IV. PHOTON STATISTICS OF THE MAZER FIELD 



The steady state distribution of photons obeys the equation 

PK,«2) = (27) 

In the Umit 52 — >■ 0, the two-photon transition probabihty in Eq. (1l8[) tends to zero and therefore we can neglect 
the third and fourth terms containing Gb2,ni,n2 in the master Eq. (pq). In this case, the upper transitions a ^ bi 
behave hke two-level atoms interacting with the mode 1 of the cavity. The lower transitions bi 62 and hence the 
two-photon transitions a —> bi —> b2 are forbidden in the interaction. The steady state solution of the Eq. (l27| ) 
can then be obtained in analytical form by using the principle of detailed balance as discussed by Meyer et al p7[ . 
The photon statistics in this two-level problem is a mixture of thermal and shifted thermal distributions. 
When g2 ^ 0, both the one- and two-photon effects of atomic transitions contribute in building up the cavity field. 
The two-photon terms (third and fourth terms) in the master equation have no counterpart in the decay terms and 
therefore the equation is not solvable analytically for the steady state distribution by the principle of detailed 
balance adopted in all the previous works on micromasers. In this general case, we integrate the master equation (|2^ ) 
numerically using fourth order Runge kutta method to get the steady state solution. We do not use any decorrelation 
approximation. The photon distribution Pi(n) and P2(n) in the cavity modes 1 and 2 are obtained respectively using 

00 00 
Pi(n) =^P(n,0, P2(n) = ^P(Z,n) . (28) 
1=0 1=0 

The normalized variances of photon distribution in the two cavity modes are defined by 

= ^""^r f , a = l,2. (29) 

In Fig. 4, we present the numerical results of the photon distribution in steady state for ultra-cold, incident atoms 
{k/K << 1) by assuming equal parameters for the decay terms Ci = C2 = C, n^j = n^^ = rib- The graph shows that for 
<?i < 92, the photon statistics in mode 1 is super- Poissonian {af > 1) while that of mode 2 is sub-Poissonian (ct| < !)• 
The photon distribution for gi > 52 is identical to that of gi < 92 except that modes 1 and 2 are interchanged. When 
gi = 92, each mode of the cavity field exhibits Poissonian statistics (cr^ « 1) of mean r/2C — 1. To understand these 



numerical results, we now approximate the master equation ( |26| ) by dropping the first and second terms corresponding 
to one-photon transitions. In fact, for the parameters of Fig. 4, the barrier and well amplitudes are p^_^ „^ « —1, 
w for wide range of rii and 712 values. Therefore, the one-photon transition probability Pm^ngC'^ ^ ^1) 
in Eq. ( [l7| ) is approximately zero. The two-photon emission probability P^i, 713(0 ^ ^2) in Eq. ( p^ ) can then be 
approximated to be 2gfg2{ni + l){n2 + l)/{gi{ni + 1) -I- 5|(ri2 + 1))^- Note that this approximation is also consistent 
with the results for ultra-cold atoms in Fig. 3. With these substitutions for the photon emission probabilities in 
gain coefficients, numerical integration of the master Eq. (^6|) again gives the same results. Moreover, in the absence 
of one-photon terms, the master equation is symmetric with respect to the labels 1 and 2 of the two cavity modes. 
Thus, the Poissonian distribution of photons in each cavity mode is purely the effect of two-photon transitions of the 
excited atoms. It should be emphasized that the transmission of atoms occurs only in the dark eigenstate component 
of the interaction during the field buildup in the cavity. The atoms interacting with the barrier-well component of 
the dressed states get reflected always. But both the reflected and the transmitted atoms have equal probability 
of two-photon emissions into the cavity. Since the one-photon emissions are forbidden in the interaction for ultracold 
atoms, the photon distribution of both the cavity modes peaks around the same photon number in Fig. 4. In general, 
the photon distribution for mode 1 peaks at a higher photon number when compared with that of mode 2 . This is 
because the mode 1 of the cavity can be populated by both one- and two-photon emissions while the mode 2 can be 
populated by only two-photon emissions from the incident atoms. We have found this behavior of steady state photon 
distribution in the case of fast, incident atoms. In Fig. 5, we display the steady state photon distribution when the 
micromaser is pumped by fast, incident atoms [k/n » 1) for the same parameters of Fig. 4(a). The graph shows 
that the fleld 1 is amplifled more than the field 2 by the stimulated, photon emissions from the incident atoms. It 
is important to note that the field 2 has an infiuence on the field 1 in the cavity through the interaction with the 
atoms. The effects of field 2 such as gain enhancement or gain reduction on field 1 have been already discussed in 
a different context viz in a two-beam laser operating on cascade three-level atoms |2^. Next, we show the effect of 
thermal photons in the cavity on the steady state photon distribution in Fig. 6. For comparison, we have also plotted 
the photon distribution for the two-level problem (^2 = 0). The length kL of the cavity is chosen to be at a resonance 
of the one-photon emission probability for the initial, excited state of the two-level atom and three photons in mode 
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1 of the cavity. The photon distribution in the two- level problem then looks similar to a mixture of thermal and 
shifted thermal distributions as discussed by Meyer et al For three-level atoms, comparison with Fig. 4(a) shows 
that the photon distribution broadens because of the presence of thermal photons {nt ^ 0) even though qualitative 
features are very similar. In particular, the two-photon emissions from the pumping atoms are still the dominant 
contribution to the steady state field for the chosen length of the cavity. The Poissonian-like statistics of photons 
in the micromaser cavity pumped by cold atoms, resembles closely the behavior of a laser operating at far above 
threshold Finally, we note that the two-photon effects are dominant over the one-photon transitions only in the 
limit (fc/K << 1) of ultra-slow motion of the incident atoms. The competition of the one- photon with the two-photon 
processes becomes stronger even for energies of the incident atoms (fc/w « 1) close to the vacuum coupling energy. 
In Fig. 7, we plot the photon distribution in the cavity for the mean momentum fc/re = 1.1 of the incident atoms. 
The inset of the Fig. 7 shows the photon emission probabilities of an excited atom for an initial, vacuum state of the 
cavity field. The graph shows that the one-photon effects participate in the field build up of the cavity and these lead 
to unequal, average number of photons in the two modes of the mazer field. 



V. SUMMARY 



We considered the interaction of mono-energetic beam of ultra-cold, cascade three- level atoms in the excited state 
with a bimodal cavity. The atom-field interaction is equivalent to reflection and transmission of the atoms through 
the potentials induced by the dressed states. There is also a reflection-less transmission of the atoms in the dark 
eigenstate of atom-field interaction. We find that two-photon transition can always be induced in an excited atom 
when its one-photon transition is forbidden. In general, the two-photon emissions from the excited atoms dominate 
over the one-photon emissions in building up the steady state field of the cavity. The steady state photon distribution 
in each mode exhibits sub- and super- Poissonian behaviors depending on the strengths of the atom-field couplings. 
The photon distribution approaches the Poissonian statistics when the atom-field coupling strengths of the two modes 
are equal. We have also obtained the photon distribution in steady state when the micromaser is pumped by fast, 
moving atoms instead of ultracold atoms. In this case, both one- and two- photon emissions from the incident atoms 
contribute in building up the steady field of the cavity. 
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FIG. 1. The scheme of the two-mode micromaser cavity pumped by S-type three-level atoms in the excited state. 
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FIG. 2. Schematic representation of the energy E of the excited atoms incident upon a two- mode micromaser 
cavity with (n, m) photons. The atom-field interaction creates a barrier (dashed) and well (dotted) potentials with 
a potential energy V = h\/g'f{n + 1) + g'H'n^ + 1) in the dressed states m+i)- '^^^ scattering from these cavity 
induced potentials leads to reflection or transmission of the atoms through the cavity. The interaction also induces a 
reflection-less transmission of the atoms in the dark state \4>n+i m+i)- However, the reflection or transmission of the 
atoms can occur only in either of the three states \a, n, m), |6i, n -|- 1, m), and \b2,n -|- 1, m -|- 1). 
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FIG. 3. The probabilities of a — > 6i (solid curve) and a ^ &i ^ &2 (dotted curve) transitions of an excited atom as a 
function of the length kL of the cavity. The cavity is initially in vacuum state and the parameters used are 92/91 = 2, 
k/K = 0.01 [(a)], k/n — 100 [(b)]. The dashed curve in Fig. 3(a) represents the photon emission probability of an 
excited two-level atom resonant with the upper transition when 92 = 0. Actual values of the dashed curve are 2.5 
times those shown. For clarity, the dotted curve in Fig. 3(a) has been displaced by 0.1 units along the Y axis. 
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FIG. 4. The steady state distribution of photons in mode 1 (sohd curve) and mode 2 (dashed curve). The parameters 
used are Ci = C2 = C, = rib^ = rib, r/C = 50, Ub = 0, kL = 200007r, k/K = 0.01 and (a) 52/51 = 2, (b) 52/51 = 1- 
In the case of 52 = 5i, the dashed curve is hardly distinguishable from the solid curve. The photon statistics for the 
parameter 52/51 = 0.5 is approximately similar to that of 52/51 = 2 except that the solid (dashed) curve corresponds 
to the photon distribution in mode 2 (1). 
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FIG. 5. The steady state distribution of photons in mode 1 (sohd) and mode 2 (dashed) for the same parameters 
of Fig. 4(a) with k/K= 100. 




FIG. 6. The steady state distribution of photons in mode 1 (sohd curve) and mode 2 (dashed curve). The parameters 
for the calculation are Ci = C2 = C, = rife^ = n^, r/C = 50, k/K = 0.01, 92/91 = 2, kL = 400007r/ and m = 1. 
The dotted curve represents the photon distribution in mode 1 for the two-level problem when 92 = 0. Actual values 
of the dotted curve are 5 times those shown. 
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FIG. 7. The steady state distribution of photons in mode 1 (solid) and mode 2 (dashed) for the same parameters 
of Fig. 4(a) with k/n = 1.1. The sohd (dashed) curves in the inset represents the one-photon (two-photon) transition 
probabiUties of an excited atom for an initial vacuum field in the cavity. The parameters used for the inset are 
52/51 = 2 and k/n = 1.1. 
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